Quantum coherence of Hard-Core-Bosons and Fermions : Extended, Glassy and Mott Phases 



O 

o 

(N 

c 

a 



Ana Maria Rey 1 - 2 , Indubala I. Satija 2 3 and Charles W. Clark 2 

1 Institute for Theoretical Atomic, Molecular and Optical Physics, Cambridge, MA, 02138, USA 
2 National Institute of Standards and Technology, Gaithersburg MD, 20899, USA and 
3 Dept. of Phys., George Mason U, Fairfax, VA, 22030.USA 
(Dated: February 6, 2008) 

We use Hanbury-Brown-Twiss interferometry (HBTI) to study various quantum phases of hard core bosons 
(HCBs) and ideal fermions confined in a one-dimensional quasi-periodic (QP) potential. For HCBs, the QP 
potential induces a cascade of Mott-like band-insulator phases in the extended regime, in addition to the Mott 
insulator, Bose glass, and superfluid phases. At critical filling factors, the appearance of these insulating phases 
is heralded by a peak to dip transition in the interferogram, which reflects the fermionic aspect of HCBs. On 
the other hand, ideal fermions in the extended phase display various complexities of incommensurate structures 
such as devil's staircases and Arnold tongues. In the localized phase, the HCB and the fermion correlations are 
identical except for the sign of the peaks. Finally, we demonstrate that HBTI provides an effective method to 
distinguish Mott and glassy phases. 
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Recent developments in the manipulations of ultra-cold 
atoms in optical lattices have opened new possibilities for 
exploring the richness and complexity of strongly-correlated 
systems. A particular topic of continuous theoretical interest 
in condensed matter physics has been the different phases in- 
duced by the interplay between disorder and interactions(lJI2] • 
In cold atomic gases disorder can be introduced in a controlled 
manner either by placing a speckle pattern on the main lattice 
or by imposing a quasiperiodic (QP) sinusoidal modulation on 
the original lattice by using an additional weaker lattice with 
the desired wavelength ratio |3]. Furthermore, by changing 
the depth of the optical potential it is experimentally possi- 
ble to control the effective dimensionality as well as the ratio 
between atomic kinetic and interaction energy in these sys- 
tems. All these attributes have made cold atomic systems a 
unique new laboratory to test many established results, ex- 
plore new phenomena underlying disordered systems and to 
confront open questions |4j|5|]. 

In this paper, we use two and four point correlation func- 
tions to compare and contrast the different quantum phases in 
a gas of one dimensional(lD) hard core bosons (HCBs) and 
fermions subject to a QP potential. In an HCB gas Ja], the 
repulsive interactions between bosons mimic the Pauli exclu- 
sion principle, and as a result HCBs resemble in many respects 
a system of non-interacting fermions. A ID HCB gas in a 
periodic optical lattice has just recently been experimentally 
realized |7]. A system of non-interacting fermions might be 
realized experimentally by reducing interatomic interactions 
via Feshbach resonances |8j. Here we treat the case of atoms 
confined by a periodic lattice, with an additional QP poten- 
tial introduced to add pseudo-random disorder. A QP system 
is in between periodic and random systems and is known to 
exhibit localization transitions at a finite depth of the addi- 
tional lattice 1 9]. Here we show that the interplay between the 
effects of disorder, interactions and quantum statistics leads 
to new quantum phases, fractal structures such as the devil's 
staircase, and fragmented Fermi distribution functions. Our 
studies reveal these effects are manifested in first- and second- 
order interferometry, which makes them accessible to direct 



experimental observation. Furthermore, we find that QP dis- 
order can be exploited to distinguish Mott from glassy phases. 

Four-point correlations can be experimentally probed by 
Hanbury-Brown-Twiss interferometry (HBTI) 1 10] which 
measures the intrinsic quantum noise in intensity measure- 
ments, commonly refereed to as shot noise or noise corre- 
lations. HBTI is emerging as one of the most important 
tools to provide information beyond that offered by standard 
momentum-distribution based characterization of phase co- 
herence. The latter is obtained from images of the density 
distribution of the atomic sample after its release from the con- 
fining potential. HBTI is performed by measuring the density- 
density correlations in such images fill H3l . Sharp peaks 
(dips) in these correlations reflect bunching (anti-bunching) 
of bosons (fermions), characterize the underlying statistics, 
and reveal information about the spatial order in the lattice 
00. 

In our study we find that the QP potential induces in the 
HBTI pattern an hierarchical set of peaks appearing at the re- 
ciprocal lattice vectors of both lattices with competing period- 
icities. We use these peaks together with the first order Bragg 
peaks in the momentum distribution to provide a detailed char- 
acterization of the different many-body phases. Three charac- 
teristic phases are known in disordered HCB systems 0: 
i) the incompressible Mott Insulator (MI) phase, ii) the super- 
fluid phase, and iii) the insulating but compressible Bose-glass 
(BG) phase. Here we also find a number of Mott-like band 
insulator phases. These phases correspond to the filling of 
the sub-bands of a fragmented energy spectrum that emerges 
from the QP superfluid phase. Transitions to these phases are 
signaled by a decrease in the intensities of the first and second 
order coherent peaks. They are found to be followed by a peak 
to dip transition in the HBTI pattern when an empty band is 
populated with few atoms. 

Ideal fermions in the localized MI and glassy phases exhibit 
first- and second-order interference patterns similar to the 
ones of HCBs, except for a different sign of the peaks. How- 
ever, in the extended phase, fermions and HCBs behave quite 
differently, with fermion interference patterns displaying vari- 
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ous complexities of incommensurate structures. When plotted 
as a function of the filling factor, their quasi-momentum distri- 
bution displays an Arnold tongue-like structure; the intensity 
of HBTI peaks exhibits a step-like pattern which resembles a 
devil's staircase 1 1511. Thus, our study points out the poten- 
tial of ultra-cold gases to provide laboratory demonstrations 
of non-linear and multi-fractal phenomena. 

In a typical experiment, atoms are released by turning off 
the external potentials. The atomic cloud expands, and is pho- 
tographed after it enters the ballistic regime. Assuming that 
the atoms are noninteracting from the time of release, prop- 
erties of the initial state can be inferred from the spatial im- 
ages 1 1 ll 1 1 3l 1 1 611 : the column density distribution image re- 
flects the initial quasi-momentum distribution, n(Q), and the 
density fluctuations, namely the noise correlations, reflect the 
quasi-momentum fluctuations, A(Q. Q'), 



h{Q) = I^e^-^}*, 



(1) 



A(Q,Q') = (n(Q)h(Q'))-(h(Q))(h(Q')). (2) 

In Eq. (0 we have assumed that both, Q, Q 1 lie inside the 
first Brillouin zone. In this paper we focus on the quantity 
A(Q, 0) = A(Q). L is the number of lattice sites and a the 
lattice constant. is a bosonic or fermionic annihilator op- 
erator at the site j. The low energy physics of a ID gas of 
strongly correlated bosons in an optical lattice modulated by 
a QP potential is well described by the Hamiltonian: 



F<™> = -j£(8J&i +1 + SJ +1 Si) + £*S' 



(3) 



where bj is the annihilation operator at the lattice site j; it 
satisfies bosonic commutation relations, plus the on-site con- 
dition = b^ 2 = which suppresses multiple occupancies. 
Here J is the hopping energy between adjacent sites, and we 
have introduced Vj — 2A cos(27rer/) to describe the additional 
QP potential. The parameter A is proportional to the intensity 
of the lasers used to create the QP lattice |3[. Note that the 
primary optical lattice defines the tight-binding condition, and 
the QP potential is created by a secondary optical lattice. Ideal 
fermions are described by the same Hamiltonian, in which the 
HCB operators are replaced by fermion operators. 

For a single atom, the eigenfunctions at site j, ip\ n ' and 
eigenenergies E^v of the Hamiltonian in Eq.0 satisfy: 

■2Xcos(2Traj)^ ( j n) = E {n) ipf ] . (4) 



This is the Harper equation, a paradigm in 1 D QP systems 1 9 ] . 
For irrational a, the model exhibits a transition from extended 
to localized states at A c = J. Below criticality, all states are 
extended Bloch-like states characteristic of a periodic poten- 
tial. Above criticality, all states are exponentially localized 
and the spectrum is point-like. At criticality the spectrum is a 
Cantor set, and the gaps form a devil's staircase. 



In this letter we treat the case ofer = (V5-l)/2. In numer- 
ical studies, a is approximated by the ratio of two Fibonacci 
numbers F M -i/F M , (F 1 = F = l,F i+1 = F % + -Fj-i), 
which describe the best rational approximant to a. For this 
rational approximation the unit cell has length Fm and the 
single-particle spectrum consists of Fm bands and Fm — 1 
gaps. The gaps occur at half-values of reciprocal lattice vec- 
tors, which we label by the Bloch index Q n = ±ncr(mod 
l)27r/a, with n an integer and Q n a 6 (— 7r, tt). The size of 
the gaps decreases as n increases. In this paper, we present 
our results for a — 55/89. 

In general the phase diagram of interacting bosons with dis- 
order is rather complicated. However the analysis simplifies 
in the HCB limit, because multiple site occupancies are for- 
bidden and the system is thus "fermionized". The ground 
state energy of an iV-particle HCB system is the sum of the 
first N single-particle eigenstates, as is the case for ideal 
fermions. This is why the localization transition occurs at 
A c = J in HCB systems. Two-point and four-point corre- 
lation functions for HCBs can be expressed as a Toplitz-like 
determinant involving two-point free fermionic propagators: 
gim = J2n=o ih Ym ^ ■ However, in contrast to the ideal 
fermion case where a direct application of Wick's theorem 
can be used to express any correlation function in terms of 
gim, for HCBs the calculations are more elaborate because 
they require the evaluation of multiple determinants I13U17I1 . 




0.1 








0.08 


| v=0.3 






0.06 


1 v=i 






0.04 


■ ■ analytic 






0.02 

















-0.02 













FIG. 1: (color online) Correlations of HCBs in the glassy phase (A = 
2 J) Top: Momentum distribution. Bottom: noise correlations. The 
central peak at Q — is truncated to highlight the QP peaks. 

We now discuss the phase diagrams for HCB and fermion 
systemes, starting with the localized phase. In the limit A 3> 
A c , the single particle wave functions are localized at individ- 
ual lattice sites, -> 6{j-l n ) and E^ -> 2A cos(2vrcr/ n ), 
where the localization center l n is such that cos(27rcr/ Il _i) < 
cos(27rcr/ n ). In this limit the determinants involved in the 
evaluation of HCB correlations become trivial and an analytic 
description is possible. The momentum distribution interfer- 
ence pattern flattens out for both HCBs and fermions, n(Q) = 
v with v = N/L the filling factor. The noise correlations 
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simplify to A(Q) - v8 Qfi ~ +L^L\£ t e * 

where s — for HCBs and s = 1 for fermions. Approximat- 
ing the sum by an integral, it is possible to show that A(Q) 
is nonzero only at the reciprocal lattice vectors, Q n , of the 
combined superlattice. At these points, 



9ll\ 



A(Q n ) = A n PS vS n fi- 



2v(l-s) (-l) s sin 2 [7 



L 



(Trn) 2 



(5) 



Eq. Q illuminates various important aspects of HBTI in the 
glassy phase: i) the explicit dependence of noise correlations 
on quantum statistics, since the boson peaks are positive and 
the corresponding fermion peaks are negative (except for the 
autocorrelation peak Ao); ii) the potential of HBTI as an ex- 
perimental spectroscopy tool, since the intensity of the n = 1 
peak is directly related to the ground state energy of the many- 
body system, Ai = ^2^=0 cos(2nnalk) — E/(L2X); iii) 
the capability of noise-correlations to clearly distinguish be- 
tween BG and MI phases, since in the Mott phase, v = 1, 
the whole hierarchical set of QP peaks disappears; iv) the uni- 
versal properties of the peaks whose intensities approach an 
asymptotic value independently of the underlying commuta- 
tion relations (except for the A peak which is always larger 
for HCBs). This is a unique feature of the BG phase, as in 
the extended phase fermionic dips are weaker than the corre- 
sponding bosonic peaks (see below). 

Fig.l shows the momentum distribution, n(Q) and noise 
correlations, A(Q), in the localized phase. Both n(Q) and 
A(Q) were computed numerically for finite A using the the- 
oretical framework discussed in Refs. fl3i [nil . It should 
be noted that even though the analytic results were derived 
for A — > oo limit, they provide a fair description of the nu- 
merical results for finite A. The main difference between the 
A = 2,1 HBTI pattern and the analytic solution is the negative 
background which is larger at low quasi-momenta. The back- 
ground reflects the fact that the single particle localized wave 
functions have finite localization length £ = l/log(A). This 
also gives rise to the finite width in the momentum distribu- 
tion, which is proportional to £. As shown in the figure, QP- 
induced localization may lead to Friedel oscillations |18] in 
the momentum distribution. These oscillations are due to the 
two-fold degeneracy of the eigenstates and can be understood 
in the strongly disordered limit by taking the eigenfunctions 
to be a superposition of the two degenerate wave functions: 
one localized at If. and the other at L — If., For an odd number 
of atoms, this leads to n(Q) = 1 + 4 cos(2Qal]y-i). 

In the extended phase two-point fermion correlations gi m 
are long-range, the Toplitz-like determinants are complicated, 
and we have found an analytical analysis to be difficult. In 
this phase our numerically obtained results are summarized 
in Figs. 2 and 3. In contrast to the glassy phase, the mo- 
mentum distribution exhibits interference peaks which signal 
the quasi-long-range coherence characteristic of this phase: a 
large peak at Q = and quasiperiodicity induced peaks at 
reciprocal lattice vectors ±hQ n whose intensity depends on 
the filling factor. These peaks also exist in the noise inter- 
ference pattern where they are narrower and are accompanied 




FIG. 2: (color online)HCBs in the superfluid phase, A = 0.5J, top: 
Momentum distribution and bottom: noise correlations where the in- 
tensities of the central peaks are scaled by a factor of 1/10 to show 
the satellite dips. 




FIG. 3: (color online) HCBs noise correlation for A = 0.5J. In 
the inset we plot the visibility for the dominant peak (see text) as a 
function of v. 



by satellite dips, immersed in a negative background. The 
later reflects the long range coherence in the extended phase 
as discussed in Ref.| 13]. It should be noted that, in contrast to 
glassy phase, few QP related peaks are visible. 

A striking aspect of this phase is the cascade of lobes (seen 
in the Fig. 2) describing a series of transitions to insulat- 
ing phases. These transitions occur at filling factors v£ 1 ^ = 
?ier(mod 1) and v c = 1 — rt(r(mod 1), which correspond to 
the fillings at which the various sub-bands underlying the QP 
spectrum are completely filled. At these critical fillings, both 
first and second order correlations depict a reduction in the 
central peak and dip intensities and also the disappearance of 
the QP induced peaks as shown in Fig. 2. Explicit computa- 
tion of the compressibility shows that the system is in fact in- 
compressible at these critical fillings. These features are rem- 
iniscent of the Mott insulating phase at v = 1, so we refer to 
these phases as Mott-like band-insulator phases. However, in 
contrast to the Mott phase, number fluctuations in these phases 
do not completely vanish, but are only somewhat reduced. 

A rather interesting consequence of the Mott-like band- 
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insulating transitions is the fact that as v is increased be- 
yond the critical value, the peaks at the reciprocal lattice vec- 
tors associated with the corresponding filled band become 
dips. This change in the sign of the peaks beyond the crit- 
ical filling can be seen by plotting the fringe visibility d n — 
A(Q„)~[A(Q n +6Q)-A(Q n -6Q)}/2, with<5Q = in/aL. 
Fig. 3 illustrates this for the dominant peaks that occur at 
±Q\. Here the peak to dip transition at v^, is signaled by a 
jump in the visibility from positive to negative. This dip in the 
second order interference pattern can be interpreted as a man- 
ifestation of fermionization. Physical insight into this behav- 
ior can be gained by looking at the noise correlation pattern 
of a single particle which can be shown to exhibit negative 
fermion-type dips: A(Q) = |zo| 2( 5qo — |^qzo| 2 with Z Q is 
the Fourier-transform of the ground state wave function. For 
HCBs these negative fringes survive even for two atoms. This 
implies that the peak to dip transition can be understood as the 
consequence of an empty band populated with few atoms. 

In contrast to the superfluid behavior of HCBs, fermions in 
the extended phase have metallic properties and different in- 
terferometric pattern. For A = 0, the sharply peaked HCB 
momentum distribution is replaced by a step-function profile 
for fermions: n(Q) = 1 for \Q\ < Qp and n{Q) = for 
101 > Qf, with Qp = nv/a being the Fermi momentum. 
For < A < A c this distribution retains part of the step- 
like profile of the free fermion gas. However, the Fermi sea 
gets fragmented and additional step-function structures cen- 
tered at different reciprocal lattice vectors of the second lattice 
develop (See Fig. 4). We call the filled states centered around 
Q = the main Fermi sea and those around the QP related 
reciprocal lattice vectors the quasi-Fermi seas. 

The quasi-Fermi seas emanating from the QP reciprocal 
lattice vectors form tongue-like structures which resemble in 
many respects the Arnold tongues that describe mode-locked 
periodic windows in the circle map As the number 

of atoms increases, the width of the tongues increase. Pre- 
cisely at the critical fillings fc and v[ n \ at which a quasi- 
band is filled up and the system becomes a band-insulator, a 
quasi-Fermi sea merges with the main Fermi sea. As A in- 
creases, more tongues become visible, and with the increase 
of v they overlap in a complex pattern. As seen in Fig.4 , we 
have tongues corresponding to both the particles and holes, 
which together describe the fragmented Fermi-Dirac distribu- 
tion when both Q and v vary. 

Noise correlations for fermions ( Fig. 4) also bear a strik- 
ing contrast to those of HCBs ( Fig. 2 and 3). One sees a 
series of plateaus at reciprocal lattice vectors, as v is varied. 
The origin of this step-like structure, where jumps occur at the 
filling factors — \Q n \a/TT, can be understood from a per- 
turbative argument as follows. For non-interacting fermions, 

for Q f 0, A(Q,0) = - ^=0 I^W" with Z Q ] bein 8 
the Fourier transform of the k th single-particle eigenfunction. 
For A = , this overlap between any two different Fourier 
components is always zero as only the ground (k — 0) state 
has a zero quasi-momentum component, ie, Zq ^ = 1. For 



small A << J, first order perturbation theory explains the 
single step observed at v$ = \Qi\a/-r:, as only Zq with 

k\ = Uj^L is nonzero. As A increases further, more steps 
are observed as other eigenvectors, acquire a nonzero Q = 
component. The steps can also be seen in the Q = plane of 
n{Q). At criticality, the whole hierarchy of steps resembles a 
devil's staircase, mimicking the multi-fractal structure of the 
gaps in an energy spectrum that is a Cantor set 0. 




FIG. 4: (color online) Fermions in the metallic phase, A = 0.5J. 
Top: Momentum distribution and Bottom: noise correlations. 

In summary, quasiperiodicity can be exploited to distin- 
guish glassy and Mott phases and induce a series of Mott-like 
band insulating phases. Furthermore, it provides an explicit 
probe to witness fermionization of bosons, and facilitates the 
experimental observation of complex fractal structures. 
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